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Jj 要 : 基于 CASA 模 型 和 Miami 模 型 计算 中 亚 五 国 (哈萨克 斯 坦 .塔吉克 斯 坦 ,吉尔吉斯 斯 坦 .土库曼 斯 坦 和 乌 
效 别 克 斯 坦 ) 的 草地 净 初 级 生产 力 (NPP ), 旨 在 探究 2001 一 2019 年 中 亚 草地 NPP 的 时 空 动态 及 其 驱动 因素 。 结 果 表 
明 :2001 一 2019 年 中 亚 地 区 草地 的 实际 净 初 级 生产 力 (ANPP) 均 值 为 118.9 g Cm?*…a'; 在 空间 尺度 上 ,平原 地 区 的 
草地 ANPP 随 纬度 增加 而 呈现 由 南 到 北 逐 渐 增 加 的 趋势 ,从 平原 到 山地 草地 ANPP 随 海拔 升 高 呈现 明显 的 增加 趋 
势 ; 在 时 间 尺 度 上 ,ANPP 总 体 旦 现 不 显著 的 增长 趋势 , 近 60% 区 域 的 草地 ANPP 未 来 的 变化 趋势 将 与 过 去 相反 ,但 
不 显著 ;在 中 亚 草地 ANPP 增 加 的 区 域 , 人 类 活动 占 主导 地 位 ;在 中 亚 草地 ANPP 降 低 的 区 域 ,气候 变化 是 主要 驱动 
力 ; 中 亚 草 地 的 ANPP 与 降水 呈现 显著 正 相关 的 区 域 达到 67.8% ,温度 与 中 亚 草地 ANPP 的 变化 关系 不 显著 ,降水 是 


影响 中 亚 草地 ANPP 变 化 的 主要 气候 因子 。 
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亚洲 中 部 干旱 区 是 全 球 最 大 的 非 地带 性 干旱 
区 之 一 ,中 亚 五 国 (哈萨克 斯 坦 吉尔吉斯 斯 坦 LER 
吉 克 斯 坦 .土库曼 斯 坦 . 乌 效 别克 斯 坦 ) 是 其 中 十 分 
重要 的 组 成 部 分 ""。 中 亚 地 区 地 表 植 被 稀 足 ,水 资 
源 严重 匮乏 ,全 球 变 暧 背景 下 ,气候 变化 剧烈 , 生 
态 系统 脆弱 5。20 世 纪 80 年 代 至 今 , 中 亚 地 区 气温 
显著 上 升 ,降水 的 地 域 性 差异 也 在 变 强 上 5 。 在 过 去 
半 个 世纪 ,气候 变化 剧烈 ,甚至 出 现 极端 气候 ,包括 中 
亚 地 区 在 内 的 全 球 各 地 的 干旱 现象 在 不 断 增 加 *。 
近 30 a, 中 亚 部 分 地 区 生态 环境 退化 ,干旱 、 草 原 退 
化 和 湖泊 萎缩 等 现象 加 剧 ,使 中 亚 地 区 成 为 亚洲 万 
至 世界 上 生态 恶化 最 严重 的 区 域 之 一 ”"。 特 别 是 21 
世纪 以 来 ,中 亚 地 区 生态 环境 恶化 较为 严重 ,在 国 
际 社会 上 引起 高 度 关注 1。 

草地 是 一 种 可 再 生 资源 ,在 全 球 陆地 植被 中 占 
有 较 大 的 面积 ,具有 维持 全 球 生态 平衡 的 重要 价 
值 ,也 是 畜牧 业 的 生产 资料 和 经 营 对 象 ,同时 ,草地 
作为 陆地 生态 系统 中 最 重要 的 组 成 部 分 之 一 ,对 全 
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球 的 碳 平衡 具有 重要 的 作用 “，。 和 气候 变化 和 人 类 
活动 的 不 利 影响 ,将 会 导致 草地 生态 系统 环境 亚 
化 ,进而 导致 草地 退化 汪 。 草 地 在 中 亚 土 地 利用 类 
型 中 占有 绝对 优势 " ,在 全 球 碳 循环 .区 域 牛 态 安 全 
和 社会 稳定 中 发 挥 着 关键 作用 。 

植被 净 初 级 生产 力 (NPP) 是 指 在 单位 时 间 和 单 
位 面积 内 ,绿色 植物 所 积累 的 有 机 质 和 干 物质 总 
量 , 对 应 了 植被 的 自然 生产 能 力 ,在 调节 全 球 碳 循 
环 与 生态 过 程 环节 中 发 挥 了 重要 作用 。 草 地 生 
态 系统 的 净 初 级 生产 力 及 其 变化 是 衡量 对 气候 因 
子 高 度 敏 感 的 生态 系统 功能 的 重要 指标 , 同时 草地 
NPP 对 非 气 候 因 素 的 变化 也 会 做 出 响应 "1。 

脆弱 的 中 亚 草 地 生态 系统 容易 受到 人 类 活动 
和 全 球 气 候 变化 的 调控 。 过 去 几 十 年 ,气候 变化 和 
人 类 活动 对 中 亚 草 地 的 影响 在 加 剧 。 研 究 表明 ,中 
地 区 的 气温 以 0.4 C (10a)"' 的 速度 在 上 升 , 随 之 
而 来 的 就 是 干旱 程度 加 剧 王 ,从 而 加 速 草地 退化 和 
荒漠 化 过 程 乓 。 同 时 ,中 亚 的 体制 和 社会 经 济 活动 
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也 发 生 了 深刻 的 变化 ,伴随 1991 年 苏联 的 解体 ”， 
该 地 区 的 草原 生态 系统 也 更 容易 受到 不 同 程度 的 
人 为 干扰 。 例 如 ,放牧 活动 的 加 剧 导 致 草原 严重 退 
化 ,大 量 奔 耕地 被 改造 成 草原 。 在 这 种 背景 下 ,中 
亚 已 经 成 为 理解 气候 和 人 类 因素 对 草原 生态 系统 
影响 的 热点 地 区 。 

在 "一带 一 路 "倡议 的 大 背景 下 ,厘清 气候 变化 
和 人 类 活动 对 中 亚 草 地 的 作用 成 了 或 需 解 决 的 问 
题 。 在 以 往 的 研究 中 ,缺少 定量 的 方法 评价 中 亚 草 
地 生产 力 的 驱动 机 制 ,讨论 人 为 因素 对 草地 生产 力 
的 影响 也 较 少 。 因 此 ,需要 一 种 定量 方法 来 评估 气 
候 因子 以 及 人 类 活动 因素 对 中 亚 草地 生产 力 的 相 
对 影响 。 本 研究 模拟 中 亚 草地 的 净 初 级 生产 力 , 月 
的 是 探究 近 20 a 中 亚 草 地 NPP 变 化 格局 ,并 且 预 测 
中 亚 草地 NPP 未 来 变化 趋势 ,同时 定量 分 析 气 候 变 
化 和 人 类 活动 在 中 亚 草 地 生产 力 变化 过 程 中 的 作 
用 。 本 文 将 为 了 解 中 亚 地 区 草地 生态 系统 动态 提 
供 有 效 信 息 ,对 未 来 草地 资源 的 合理 利用 提供 
参考 。 


1 材料 与 方法 


1.1 研究 区 概况 

中 亚 地 处 欧 亚 大 陆 腹地 ,东西 横 跨 46°29’47”~ 
87°18'55”E, PA JE JES 35°07'43"~55°26'28"N , P5 2) 
里 海 和 伏尔加 河 , 东 与 中 国 接壤 , 北 到 西伯 利 亚 大 
草原 , 南 与 伊朗 、 阿 富 汗 交界 。 总 面积 约 4.0x105 
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尔 盖 台 地 和 哈萨克 斯 坦 丘 陵 ,东部 和 东南 部 是 天 山 
山脉 和 帕 米 尔 高 原 。 中 亚 地 区 海拔 自 西部 平原 向 
东部 山区 逐渐 升 高 ,表现 为 东南 高 西北 低 ,主要 地 
形 为 平原 .山地 及 丘陵 ,沙漠 面积 广 。 主 要 气候 类 
型 为 温带 大 陆 性 气候 ,包括 温带 草原 ,高原 山地 及 
温带 荒漠 气候 。 中 亚 地 区 年 均 降 水 相对 稀少 ,不 同 
地 域 的 年 均 降 水 量 相 差 大 ;日 照 充足 ,气温 日 较 差 
很 大 ,从 多 年 平均 气温 分 布 来 看 ,西部 和 南部 较 
高 。20 世 纪 80 年 代 至 今 ,中 亚 地 区 气温 显著 上 升 ， 
但 降水 存在 较 强 的 地 域 性 差异 ”。 中 亚 的 主要 植 
被 类 型 为 林地 .草地 和 耕地 ,草地 分 布 最 广泛 ,面积 
也 相对 较 大 。 

1.2 数据 来 源 

1.2.1 NPP、NDVI 及 土地 利用 数据 利用 MOD13A1 
产品 作为 NDVI 的 源 数据 ,空间 分 辨 紊 为 500 mx 
500 m, 时 间 分 辩 率 为 16 d, 采 用 最 大 化 合成 法 (max- 
imum value composite , M VC ) 将 源 数据 的 分 辩 率 合成 
为 30 d。 利 用 MOD17A3H 产品 作为 NPP 的 验证 数 
据 , 空 间 分 辩 率 为 500 mx500 mo MCD12Q1 产品 作 
为 土地 利用 数据 ,空间 分 辨 紊 为 500 mx500 m ,参考 
国际 地 圈 生 物 圈 计划 (International Geobiosphere 
Programme, IGBP) 的 土地 分 类 体系 ,将 中 亚 地 区 的 
土地 利用 类 型 重新 归 类 为 :耕地 .草地 nc LR 
PK 灌木 .建筑 用 地 .水 体 和 草地 (图 la, 其 中 草地 
包括 ICGBP 分 类 系统 中 的 典型 草原 和 稀 玻 草原 )。 
1.2.2 气象 数据 气象 数据 选用 英国 East Anglia 大 


km, FRAT AIA 57 ES 2€ , BIE Ge vo 
尔 吉 斯 斯 坦 .塔吉克 斯 坦 .土库曼 斯 坦 和 乌 效 别克 
斯 坦 。 中 亚 西部 是 图 兰 低地 ,其 北部 与 东北 部 是 图 


学 Climatic Research Unit (CRU ) 提 供 的 2001—2019 
年 的 降水 和 气温 格 点 数据 ,空间 分 辨 率 0.5°x0.5° 
(2450 km), 利 用 ArcGIS 软件 将 其 重 采样 至 500 mx 
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注 : 底 图 采用 自然 资源 部 标准 地 图 制作 , 审 图 号 为 GS(2016)1666 号 ,对 底 图 边界 无 修改 。 下 同 。 


图 1 研究 区 概况 及 采样 点 示意 图 


Fig. 1 Overview of study area and sampling point 
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500 mo CRU 通过 整合 已 有 的 若干 个 知名 数据 库 ， 
重建 了 一 套 覆 羡 完 整 的 地 表 气 候 要 素数 据 集 ,覆盖 
全 球 陆 地 的 9 个 地 表 变 量 , 具 有 较 高 可 信和 度 。 陈 发 
虎 等 2 也 认为 ,在 中 亚 地 区 1930 年 以 后 的 CRU 资料 
可 靠 性 和 适用 性 较 高 。 
1.2.3 实测 数据 所 使 用 的 实测 数据 来 自 2018 一 
2019 年 6 一 7 月 ,分别 在 哈萨克 斯 坦 .塔吉克 斯 坦 和 
吉尔 吉 斯 斯 坦 采集 的 草地 植被 生物 量 数据 。 选 取 
样 地 共 40 个 (图 1b) ,每 个 样 地 布设 5 个 1 mx] m ff 
FED ,收取 每 个 样 方 的 地 上 生物 量 (对 于 草地 中 的 
灌 丛 ,选取 当年 的 新 枝 和 叶片 作为 地 上 生物 量 ) , 求 
均值 得 到 每 个 样 地 的 生物 量 。 地 下 生物 量 参考 马 
安娜 等 所 对 中 国 草 地 植被 地 上 和 地 下 生物 量 的 关 
系 的 研究 求 得 。 
1.3 研究 方法 
1.3.1 NPP% ”本 研究 定义 了 3 种 类 型 NPP, 分 别 
是 由 CASA 模 型 计算 的 ANPP( 人 类 活动 与 气候 因素 
共同 驱动 的 NPP) ,Miami 模 型 模拟 的 气候 因素 驱动 
的 潜在 净 初 级 生产 力 (PNPP) ,以 及 人 类 活动 占用 的 
净 初 级 生产 力 (HNPP)。 
采用 CASA 模 型 来 模拟 中 亚 地 区 的 草地 ANPP 
(实际 净 初 级 生产 力 ) ,该 模型 中 植被 NPP 由 光合 有 
效 辐射 (APAR ) 和 实际 光 能 利用 率 (g) 决 定 : 
ANPP(x,1) - APAR(x,t) x &(x,t) (1) 
APAR 的 计算 公式 如 下 : 
APAR(x,t) =SOL(x,t) x FPAR(x,t)x0.5 (2) 
SUA: SOLARI A BER at ; FPAR (x,t) Ro HL 
被 层 对 入 射 光 合 有 效 辐射 的 吸收 比例 ,可 以 通过 归 
一 化 植被 指数 (NDVD) 求 得 ,常数 0.5 表 示 植 被 所 能 
利用 的 太阳 有 效 辐射 的 比例 。 
光 能 利用 率 的 计算 公式 如 下 : 
é(x,t) = T (x,t) x T (x,t) x W (x,t) X E ux (3) 
式 中 :Tax 和 7s2(x 力 表示 低温 和 高 温 的 胁迫 系数 ; 
取 ( 思 为 水 分 胁迫 系数 ;se 是 理想 条 件 下 的 最 大 光 
能 利用 率 。 
采用 Miami 模 型 模拟 中 亚 草 地 的 PNPP。Miami 
模型 根据 温度 (7) 和 降水 (P) 资 料 建立 ,广泛 用 于 计 
算 大 尺度 范围 的 气候 驱动 的 草地 的 潜在 净 初 级 生 
产 力 (PNPP) ,计算 公式 如 下 : 
PNPP = min (NPP,, NPP,) (4) 
NPP, = 3000/(1 + e? | (5) 


NPP, = 3000 x (1 — e °°") (6) 


式 中 :NPP; 和 NPP, 分 别 表示 温度 和 降水 引起 的 草地 
NPP。 
HNPP 为 PNPP 和 ANPP 的 差 值 : 
HNPP=PNPP-ANPP (7) 

1.3.2 NPP 的 验证 基于 实测 数据 开展 CASA 模型 
模拟 结果 的 验证 工作 ,需要 将 地 面 实测 生物 量 转 换 
为 NPP。 公 式 如 下 : 

NPP, = (AGB + BGB) X turnover (8) 
CUP AGB 为 地 上 生物 量 ;BCB 为 地 下 生物 量 ;turn- 
over 是 将 单位 面积 生物 量 数据 换算 为 实测 NPP 值 的 
转换 系数 ,本 研究 参考 Chen 等 "的 研究 , 取 固 定 值 
0.5。 
1.3.3 NPP 的 趋势 分 析 采用 基于 最 小 二 乘法 的 线 
性 回归 分 析 计 算 中 亚 地 区 草地 ANPP 2001—2019 4E 
年 际 变化 趋势 ,并 利用 MATLAB 软件 计算 NPP 的 线 
性 回归 系数 : 


nx Y (ixNPP)- yixyNPP 
i=l i=l izl 


nx Yi 一 ($) 

SUA : Slope 是 趋势 线 的 线性 回归 系数 , 当 Slope >0 
时 ,说 明 草 地 NPP 在 年 内 的 变化 趋势 是 增加 的 , 反 
之 则 减少 。 

1.3.4 NPP 趋 势 的 持续 性 分 析 KA Hurst 指数 分 
析 中 亚 草地 ANPP 趋势 的 持续 性 ,通过 MATLAB 软 


Slope = (9) 


件 实现 。 
均值 序列 : 
ELSEN, 1=1, 2, 3, … — (10) 
累计 离 差 : 
xar) Yle- sise a1) 
极 差 : 
R(x) 3 max X(1,z) i min X(t,7), (12) 
7:1,2,3,- 
标准 差 ; 
sof eo- as 
7-1,2,3, + 


4E X. R/S-R (rS) , £ R/S coc", HAA PAT AY EY 
间 序 列 符合 Hurst 指数 的 定义 ,五 称 为 Hurst 指数 , 范 
围 为 0~1。 若 0.5<E<1, 表 明 所 研究 的 时 间 序 列 变化 
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状态 与 过 去 一 致 , 即 正 向 持续 ; 若 0«H«0.5 , EH] Ir 
研究 时 间 序 列 的 变化 状态 与 过 去 相反 , 即 反 向 
持续 。 

1.3.5 NPP 与 温度 、 降 水 的 相关 性 分 析 利用 Pear- 
son 相关 关系 法 对 中 亚 草地 ANPP 和 中 亚 地 区 年 均 
温 与 年 降水 之 间 的 相关 关系 进行 了 分 析 ,并 通过 
MATLAB 软件 实现 。 

2 36.-7) 


R,= (14) 


n 


Ee- [6-9 


IEP in Jg] f] Fe AR BE 5 x; Aly IRR LAE I 
ANPP 值 和 第 ;年 的 年 均 温 或 年 降水 ;元 和 7 了 分 别 表 
示 ANPP 的 多 年 平均 值 和 多 年 平均 气温 或 多 年 平均 
降水 ;相关 系数 RR 的 绝对 值 的 范围 为 0~1。 


2 结果 与 分 析 


2.1 模型 检验 

利用 2018 一 2019 年 哈萨克 斯 坦 .塔吉克 斯 坦 和 
吉尔 吉 斯 斯 坦 的 生长 季 实 测 NPP, 对 CASA 模型 模 
拟 的 草地 NPP 进 行 验证 (图 2a) ,采用 的 实测 数据 样 
方 分 布 范 围 广 ,涵盖 多 种 草地 类 型 ,监测 时 间 也 较 
为 一 致 ,可 以 较 好 地 代表 中 亚 地 区 草地 净 初 级 生产 
Fi SEAM) NPP FI CASA EAA NPP PRE ABER? 
0.69 (P«0.001) ,表明 CASA 模型 是 适用 于 估算 当地 
草地 NPP。 

为 了 进一步 评估 CASA 模型 对 中 亚 地 区 NPP 的 
模拟 精度 ,利用 AreGIS 10.2 的 渔网 工具 ,在 研究 区 
内 均匀 生成 500 个 采样 点 ,提取 2001 一 2019 年 研究 
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区 内 CASA 模型 模拟 NPP 值 和 MODIS NPP 值 ,用 于 
验证 分 析 ( 图 2b)。 验 证 结果 表明 ,MODIS NPP 数 值 
略 高 于 CASA 模型 模拟 NPP, 二 者 拟 合 程 度 较 高 
(R-0.79 , P«0.001) > 
2.0 中 亚 草 地 NPP 时 空 变化 特征 及 变化 趋势 

在 海拔 和 气候 的 影响 下 ,整个 区 域 的 ANPP 具 
有 明显 的 空间 异 质 性 (图 3)。 模 型 模拟 得 到 ， 
2001 一 2019 年 ,中 亚 地 区 草地 的 年 均 ANPP 为 118.9 
gC"m**a!。 草 地 ANPP 最 高 和 最 低 的 年 份 分 别 为 
2002 年 (135.5 g C-m^?-a ^) F 2008 4E (97.6 g C-m^- 
al)。 分 析 得 出 ,中 亚 草地 的 平均 ANPP 的 低 值 区 面 
积 (0~60g Cnm”.a ) 占 中 亚 地 区 草地 面积 的 19.2%， 
主要 分 布 区 域 为 里 海 东 部 海岸 线 经 哈萨克 斯 坦 南部 
至 巴尔 喀什 湖 ; 高 值 区 面积 (>240 g Cem? a") rH PIE. 
地 区 草地 面积 的 6.0% , 主要 分 布 在 哈萨克 斯 坦 北 部 
草原 西天 山 和 帕 米 尔 高 原 西 部 。 中 亚 地 区 的 年 均 
草地 ANPP 大 部 分 集中 于 60~240 g Cm …a' 之 间 。 
从 空间 分 布 来 看 ,平原 地 区 的 草地 ANPP 随 纬度 增 
加 而 呈现 由 南 到 北 逐 渐 增 加 的 趋势 ,从 平原 到 山地 
草地 ANPP 随 海拔 升 高 呈现 明显 的 增加 趋势 。 

吉尔 吉 斯 斯 坦 和 塔吉克 斯 坦 都 是 典型 的 山地 
国家 ,草地 面积 占 两 国 国土 总 面积 的 比例 约 60%。 
2001 一 2019 年 ,两 国 的 年 均 草 地 ANPP 高 达 158.52 
gCma', 且 星 现 低 海拔 到 高 海拔 逐渐 增加 的 趋 
势 。 哈 萨克斯 坦 国土 面积 占 中 亚 五 国 总 面积 的 比 
例 超过 75%, 境 内 土地 利用 类 型 较为 丰富 ,草地 在 
土地 利用 类 型 中 占有 绝对 优势 ,该 国 的 草地 年 均 实 
际 净 初级 生产 力 为 115.93 g C-m?-a^, HÆ ANPP 
随 纬度 升 高 逐渐 增加 的 趋势 。 乌 效 别 克 斯 坦 和 十 


500 R?-0.79 e 
P<0.001 * 


模拟 NPP/(g C m? a) 


MODIS NPP/(g C : m? * a") 


图 2 CASA 模 型 模拟 NPP 与 实测 NPP 及 MODIS NPP XJ FE 
Fig.2 Comparison of CASA model simulated NPP with measured NPP and MODIS NPP 
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图 3 2001 一 2019 年 中 亚 草地 ANPP 的 平均 空间 分 布 
Fig. 3 Averaged spatial distribution of grassland ANPP in 
Central Asia from 2001 to 2019 


库 曼 斯 坦 两 国境 内 荒漠 为 主要 土地 利用 类 型 ,草地 
面积 仅 占 两 国 总 面积 的 27% ,两 国 的 年 均 草 地 
ANPP 值 相对 偏 小 , 仅 为 89.88 gC.m a'。 

回归 分 析 ( 图 4) 表 明 ,2001 一 2019 年 ANPP 增 加 
区 域 的 面积 占 中 亚 草地 面积 的 60.4% , 主要 分 布 在 
哈萨克 斯 坦 南部 和 东北 部 草原 以 及 天 山西 昔 地 区 ; 
ANPP 降 低 区 域 的 面积 占 中 亚 草 地 面积 的 39.6% , E 
要 分 布 在 哈萨克 斯 坦 西 部 和 北部 草原 。 超 过 80% 
的 草地 的 ANPP 变化 趋势 不 显著 。ANPP 趋势 通过 
显著 性 检验 (P<0.05 ) 的 区 域 面积 占 中 亚 草地 面积 的 
18.6% , 主要 集中 在 哈萨克 斯 坦 东 南部 四 州 ,其 中 
13.3% 的 区 域 草地 ANPP 显 著 增 加 ,5.3% 的 区 域 草地 
ANPP 显著 降低 。 

2001 一 2019 年 中 亚 草地 ANPP 变 化 趋势 持续 性 
分 析 (Hurst 指数 ) 的 结果 表明 (图 5),ANPP 未 来 变化 
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图 4 2001 一 2019 年 中 亚 草地 ANPP 的 变化 趋势 
Fig.4 Change trend of ANPP in Central Asia grassland 
from 2001 to 2019 
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图 5 2001 一 2019 年 中 亚 草地 ANPP 的 趋势 变化 的 持续 性 
Fig. 5 Sustainability of ANPP trend change in Central Asia 
grassland from 2001 to 2019 


趋势 与 过 去 相同 (8>0.5) 的 区 域 面 积 占 中 亚 草 地 面 
积 的 40.2% , 主要 分 布 在 哈萨克 斯 坦 ,ANPP 未 来 变 
化 趋势 与 过 去 相反 (CE<0.5) 的 区 域 面积 占 中 亚 草 地 
面积 的 59.8% , 主要 分 布 在 乌 效 别克 斯 坦 东 南部 以 
及 哈萨克 斯 坦 北 部 。 

2001—2019 年 中 亚 地 区 草地 的 年 均 PNPP 为 
466.89 g Cn”a (6a) ,中 亚 草 地 的 平均 ANPP 的 
低 值 区 面积 (<300 g C.m2.a0) 占 中 亚 地 区 草地 面积 
的 17.496 , 主要 分 布 在 哈萨克 斯 坦 的 克 孜 勒 奥 尔 达 
州 和 马 兹 别克 斯 坦 ;高 值 区 面积 (>700 g Cem- a) 
占 中 亚 草 地 面积 的 7.5% , 主要 分 布 在 塔吉克 斯 坦 和 
哈萨克 斯 坦 东 部 边境 以 及 吉尔 吉 斯 斯 坦 的 山区 。 
人 研究 发 现 , 中 亚 地 区 草地 的 PNPP 的 空间 分 布 规 律 
与 ANPP 高 度 一 致 。2001 一 2019 年 中 亚 草地 的 
PNPP 呈现 增加 趋势 的 区 域 面积 占 草 地 面积 的 
31.8% (Al 6b) , 主要 分 布 在 哈萨克 斯 坦 的 东部 和 塔 
吉 克 斯 坦 的 山区 , 除 此 之 外 ,大 部 分 区 域 的 PNPP 呈 
现 降低 趋势 , 占 草 地 面积 的 68.2%。 整 体 来 看 ,中 亚 
草地 PNPP 的 变化 趋势 并 不 显著 ,其 中 显著 降低 的 
区 域 占 比 仅 为 4.7% ,主要 分 布 在 哈萨克 斯 坦 的 西 
北部 。 

2.3 中 亚 草 地 NPP 变 化 的 驱动 因素 

将 气候 变化 和 人 类 活动 对 草地 NPP 影 响 定 义 
为 6 种 情景 ( 表 1) ,其 中 SANPP SPNPP 和 SHNPP 分 
别 代表 ANPP PNPP 和 HNPP 的 斜率 。 

2001 一 2019 年 中 亚 草 地 ANPP 增 加 的 区 域 面积 
占 比 为 60.4%( 图 7)。 其 中 由 气候 变化 导致 的 ANPP 
增加 的 区 域 面积 占 草地 面积 的 16.0% ,主要 分 布 在 
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6 2001 一 2019 年 中 亚 草 地 的 平均 PNPP 以 及 PNPP 的 变化 趋势 
Fig. 6 Average PNPP of grassland in Central Asia and its change trend from 2001 to 2019 


#1 评估 中 亚 草地 ANPP 变化 驱动 因素 的 6 种 情景 


Tab.1 Six scenarios of evaluating the driving factors of ANPP change in Central Asia grassland 


SANPP 情景 SPNPP SHNPP 定义 
情景 1(IDC) + + 气候 变化 导致 草地 ANPP 增 加 
SANPP>0 情景 2(IDH) - - 人 为 因素 导致 草地 ANPP 增 加 
情景 3(IDCH) + - 气候 变化 和 人 为 因素 共同 导致 草地 ANPP 增 加 
lise 4 (DDC) - = 气候 变化 导致 草地 ANPP 减 少 
SANPP<0 情景 5(DDH) * * 人 为 因素 导致 草地 ANPP 减 少 
情景 6(DDCH) - + 气候 变化 和 人 为 因素 共同 导致 草地 ANPP 减 少 
z OWE OB WE SE WB, 39.6%。 其 中 由 气候 变化 导致 的 ANPP 降 低 的 区 域 
3 面积 占 草地 面积 的 31.8% , 主要 分 布 在 哈萨克 斯 坦 
的 西部 和 北部 草原 ;由 人 类 活动 导致 ANPP 降 低 的 
a £ 区 域 面积 占 比 为 4.8% ,主要 分 布 在 哈萨克 斯 坦 的 西 
南部 ;由 气候 变化 和 人 类 活动 共同 导致 的 ANPP 降 
z 低 的 区 域 面积 占 比 仅 为 3.1%。 因 此 ,在 中 亚 草 地 
ai E aes S ”ANPP 增 加 的 区 域 , 人 类 活动 占 主导 地 位 ;在 中 亚 草 
情景 i 地 ANPP 降 低 的 区 域 ,气候 变化 是 主要 驱动 力 。 
mmm DC mmm DH = IDCH 0 500km 
Z [ = DDC == DDH ww DDCH 一 一 £ 为 探究 影响 2001 一 2019 年 中 亚 草地 ANPP 变 化 
T oE 60°E 70°E 80°E 90°E 的 主要 气候 因子 ,本 研究 分 别 做 了 ANPP 与 降水 (图 


图 7 气候 变化 和 人 类 活动 对 中 亚 草地 ANPP 
变化 相对 影响 的 空间 分 布 
Fig.7 Spatial distribution of relative impacts of 
climate change and human activities on grassland 
ANPP change in Central Asia 


哈萨克 斯 坦 东 部 地 区 和 帕 米 尔 高 原 的 西部 ;而 由 人 
类 活动 导致 ANPP 增 加 的 区 域 面 积 占 比 为 33.4% , E 
要 分 布 在 哈萨克 斯 坦 境内 ;气候 变化 和 人 类 活动 共 
同 导致 ANPP 增 加 的 区 域 面 积 占 比 为 11.0% ,主要 分 
布 在 哈萨克 斯 坦 东 部 地 区 和 吉尔 吉 斯 坦 。2001 一 
2019 年 中 亚 草 地 ANPP 降低 的 区 域 面 积 占 比 为 


8a) 以 及 温度 (图 8b) 的 相关 性 分 析 。2001 一 2019 
年 ,除去 天 山 山脉 和 帕 米 尔 高 原 的 部 分 地 区 ,中 亚 
绝 大 部 分 草地 的 ANPP 与 降水 呈现 正 相关 , 占 比 达 
到 96.4% ,二 者 显著 正 相关 的 区 域 面积 占 中 亚 草 地 面 
积 的 67.8%。 降 水 是 影响 中 亚 草 地 ANPP 变 化 的 主要 
气候 因子 ,温度 与 中 亚 草地 ANPP 的 变化 关系 不 
显著 。 


3 讨论 


2001 一 2019 年 ,中 亚 地 区 草地 ANPP 均 值 为 
118.9 g Cm?…a!, 本 研究 结果 小 于 Chen 等 利用 
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图 8 2001 一 2019 年 中 亚 草 地 ANPP 与 年 降水 及 年 均 温 的 相关 性 


Fig. 8 Correlation of ANPP with annual precipitation and annual mean temperature in Central Asia grassland from 2001 to 2019 


CASA 模型 模拟 的 中 亚 五 国 及 新 疆 草 地 生态 系统 的 
NPP(1982—2015 年 )。 研 究 时 段 不 同 , 人 研究 区 范围 
有 所 差异 ,采用 不 同 数 据 源 的 NDVI. 土 地 覆盖 类 型 
等 数据 ,会 造成 模拟 值 间 的 差异 。 朱 文 泉 等 3 研究 
也 表明 ,对 于 不 同 研究 区 用 于 模型 模拟 的 最 大 光 能 
利用 率 也 不 同 。 不 同 的 研究 区 以 及 不 同 的 数据 源 
所 造成 的 研究 结果 差异 ,需要 有 实测 数据 加 以 佐 
证 。 本 研究 分 别 利 用 实测 NPP 和 MODIS NPP 对 
ANPP 进 行 验证 ,二 者 与 本 研究 模拟 的 ANPP 都 有 较 
高 的 拟 合 度 。 另 外 ,CASA 模型 模拟 的 NPP(2001 一 
2019 年 ) 整 体 低 于 MODIS 的 NPP, 这 与 刘 洁 等 ”的 
研究 结论 一 致 。 

2001—2019 年 ,中 亚 大 部 分 地 区 的 草地 ANPP 
的 未 来 变化 趋势 将 与 过 去 相反 ,表明 大 部 分 地 区 的 
草地 ANPP 在 未 来 会 有 减少 趋势 ,这 可 能 是 由 于 牲 
畜 量 的 快速 增长 ,会 导致 超载 放牧 的 现象 加 剧 ,从 
而 致使 草地 生产 力 降 低 ”。 平 原 地 区 的 草地 ANPP 
随 纬度 增加 而 呈现 由 南 到 北 逐 渐 增 加 的 趋势 ,平原 
到 山地 草地 ANPP 随 海 拔 升 高 呈现 明显 的 增加 趋 
势 , 这 与 Chen 等 后 的 研究 结果 一 致 。 塔 吉 克 斯 坦 和 
吉尔 吉 斯 斯 坦 的 山区 的 草地 ANPP 明 显 高 于 中 亚 其 
他 地 区 ,山区 降水 量 高 是 其 主要 原因 。 

半 干 旱 和 干旱 生态 系统 中 的 植被 对 气候 变化 
极其 敏感 下 。2001 一 2019 年 由 气候 变化 导致 中 亚 
草地 ANPP 增 加 的 区 域 面积 占 草地 面积 的 16.0%; 由 
气候 变化 导致 草地 ANPP 降低 的 区 域 面 积 占 草地 面 
积 的 31.8%。 本 研究 发 现 , 降 水 是 影响 中 亚 草 地 
ANPP 的 主要 因素 。 相 关 人 研究 表明 ,中 亚 草 地 NPP 
对 降水 变化 比 对 气温 变化 更 敏感 加。 降水 在 很 大 
程度 上 决定 了 草地 植被 的 生长 。 降 水 减少 会 降低 


草原 植物 的 光合 作用 效率 ,抑制 植物 活动 和 有 机 质 
生产 ,最 终 降低 草地 生产 力 ; 反 过 来 ,降水 增加 可 
以 提高 半 干 旱 和 王 旱 生态 系统 中 的 土壤 水 分 含量 ， 
从 而 提高 草地 生产 力 呈 。 刘 国 华 等 后 也 认为 ,降水 
量 增加 可 以 改善 土壤 水 分 条 件 ,同时 提高 土壤 水 分 
的 有 效 性 和 土壤 水 分 利用 率 , 降 低 土壤 水 分 胁迫 ， 
从 而 有 利于 草地 植被 生长 。2001 一 2019 年 ,由 气候 
变化 导致 草地 NPP 降低 的 区 域 主要 分 布 在 哈萨克 
斯 坦 西 北部 和 南部 的 部 分 地 区 ,由 气候 变化 导致 
NPP 增 加 的 区 域 主要 分 布 在 哈萨克 斯 坦 东 部 的 部 分 
地 区 ,这 些 区 域 年 降水 量 的 变化 直接 导致 了 ANPP 
的 变化 。 另 外 ,由 于 温度 的 升 高 导致 草地 植被 气孔 
关闭 ,同时 会 伴随 着 洪 在 蒸 散 和 水 分 胁迫 增强 , 进 
而 导致 草 地 生产 力 降低 忆 , 所 以 哈萨克 斯 坦 西 部 边 
境 存在 温度 升 高 显著 影响 草地 NPP 降 低 的 现象 。 

人 类 活动 是 中 亚 草原 净 初 级 生产 力 动态 变化 
的 重要 驱动 力 呈 。2001 一 2019 年 由 人 类 活动 导致 
中 亚 草 地 ANPP 增 加 的 区 域 面 积 占 草地 面积 的 
33.4% ,主要 分 布 在 哈萨克 斯 坦 的 东部 。 这 是 由 于 苏 
联 解 体 后 经 济 体制 的 变化 ,缓解 了 过 度 放牧 ,哈萨克 
斯 坦 的 大 量 牧场 被 遗弃 并 转变 为 半 自 然 草 原 5。 同 
时 耕地 大 量 弃 耕 转变 为 草地 ,到 2010 年 , 哈 陕 克 斯 
坦 原 有 的 耕地 仍 有 58% 被 弃 耕 ,大 面积 的 农田 处 于 
闲置 状态 ,已 经 转换 为 天 然 草地 ,导致 草地 生产 力 
有 所 增加 中 。 另 外 ,2000 年 以 后 ,由 于 政策 改革 和 
政府 对 农业 的 支持 ,哈萨克 斯 坦 的 农作物 和 畜牧 业 
生产 开始 恢复 ,哈萨克 斯 坦 政府 实施 农业 多 样 化 
政策 ,以 支持 畜牧 业 生 产 ,使 得 草地 生产 规范 化 ” ， 
这 同样 也 有 利于 草地 植被 生产 力 的 增加 。Han E 
认为 干旱 地 区 合理 的 草地 放牧 制度 可 以 有 效 地 防 
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止 中 亚 草 原生 态 系统 的 退化 ,促进 草地 生产 力 的 增 
加 ,有 利于 草地 生态 系统 的 可 持续 发 展 。 相 关 研 究 
表明 ,健康 的 草地 状况 可 以 有 效 地 缓冲 放牧 和 气 
AT HC ,并 能 够 维持 较 高 的 草地 生产 力 和 生物 多 
样 性 ,而 不 放牧 并 不 一 定 能 持续 维持 草地 的 健康 
状况 3 ,适度 放牧 也 会 促进 草地 生产 力 增加 ”!。 然 
而 ,不 同 放牧 条 件 下 草原 净 初 级 生产 力 和 放牧 之 间 
的 关系 的 机 制 仍 需 展开 深入 研究 ” 。 

本 文 研究 2001—2019 年 中 亚 草地 净 初 级 生产 
力 动态 ,时 间 尺 度 相 对 较 短 ,空间 尺度 相对 较 大 。 
今后 研究 需要 获取 更 长 时 间 序 列 的 数据 ,着 重 研究 
小 尺度 区 域 , 以 求 更 加 系统 地 掌握 中 亚 草地 的 生产 
力 的 动态 变化 及 驱动 机 制 。 


4 结论 


本 研究 基于 CASA 模型 计算 2001 一 2019 年 中 亚 
地 区 草地 的 实际 净 初 级 生产 力 (ANPP) ,进而 分 析 其 
时 空 动 态 趋势 ;通过 定义 不 同 的 情景 来 分 析 人 类 活 
动 和 气候 变化 对 草地 ANPP 的 驱动 机 制 ,确定 气候 
和 人 为 因素 对 草地 ANPP 的 影响 ,得 出 以 下 主要 
结论 : 

(1) 2001 一 2019 年 中 亚 地 区 草地 的 ANPP 均 值 
为 118.9g C* m^: a” ,变化 范围 为 97.6~135.5 g C- m° 
al。 草 地 ANPP 在 空间 尺度 上 呈现 从 哈萨克 斯 坦 南 
部 到 哈萨克 斯 坦 北 部 、 从 低 海拔 地 区 到 高 海拔 地 区 
逐渐 增加 的 规律 。 

(2) 2001 一 2019 年 中 亚 地 区 草地 的 ANPP 旦 现 
增加 趋势 的 区 域 占 比 为 60.4%, 且 超过 80% 的 草地 
的 ANPP 变 化 趋势 不 显著 ,ANPP 总 体 呈 现 不 显著 的 
增长 趋势 。 接 近 60% 的 中 亚 草地 的 ANPP 未 来 会 呈 
现 不 显著 的 负 癌 持续 性 的 变化 趋势 。 

(3) 2001 一 2019 年 中 亚 地 区 草地 的 PNPP 均 值 
为 466.89 g Cem? a ',PNPP 呈 现 降 低 趋 势 的 区 域 占 
比 为 68.2% , 且 超 过 95% 的 草地 的 PNPP 变 化 趋势 不 
显著 ,PNPP 总 体 呈 现 不 显著 的 降低 趋势 。 

(4) 2001—2019 年 由 人 类 活动 导致 中 亚 草 地 
ANPP 增 加 的 区 域 面积 占 草地 面积 的 33.4%; 由 气候 
变化 导致 的 ANPP 降低 的 区 域 面 积 占 草地 面积 的 
31.8%。 在 中 亚 草地 ANPP 增 加 的 区 域 ,人 类 活动 占 
主导 地 位 。 但 是 ,在 中 亚 草地 ANPP 降 低 的 区 域 , 气 
候 变 化 是 主要 驱动 力 。 

(5) 2001—2019 年 中 亚 草 地 的 ANPP 与 降水 呈 
现 显著 正 相 关 的 区 域 占 比 达到 67.8% ,温度 与 中 亚 
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草地 ANPP 的 变化 并 无 显著 关系 。 降 水 是 影响 中 亚 
草地 ANPP 变 化 的 主要 气候 因子 。 
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Study on spatio-temporal dynamics and driving factors of 
NPP in Central Asian grassland 
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Abstract: Based on the CASA and Miami models, the net primary productivity (NPP) of five Central Asian 
countries (Kazakhstan, Tajikistan, Kyrgyzstan, Turkmenistan, and Uzbekistan) was calculated in this study. The 
purpose was to explore the Spatio-temporal dynamics and driving factors of NPP in Central Asian grassland from 
2001 to 2019. The results showed that the average ANPP of grassland in Central Asia from 2001 to 2019 was 
118.9 g C- m ^-:a '. On the spatial scale, the ANPP of grassland in the plain area increased gradually from south to 
north with increasing latitude. In contrast, the ANPP of grassland from plain to mountain grassland increased with 
increasing altitude. On the temporal scale, the overall growth trend of ANPP was not significant, and the future 
change trend of grassland ANPP in approximately 60% of the region was opposite to that of the past. The trend 
was not significant in the region of Central Asia with increased grassland ANPP where human activities are 
dominant and with decreased grassland ANPP, with climate change as the main driving force. A significant 
positive correlation existed between ANPP and precipitation accounting for 67.8% , with no significant 
relationship between temperature and the change of ANPP in Central Asian grassland. Therefore, precipitation is 
the main climatic factor affecting the change of ANPP in Central Asian grassland. 
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